COURSE-KEEPING 
WITH  AUTOMATIC  CONTROL 


Eduardo  Aguayo 


■ill  ATE 


ontepey,Gali    . 


T 


c": 


:  :.-,•_.. 


COU  RS  E  -  K  E  EPIN  G 


WITH  AUTOMATIC  CONTROL 


by 


Eduardo  Aguayo 


Thesis  Advisor: 


G.    J.    Thaler 


■BBBMHBK  ?*  .  f-.r-.:   , 


December   1973 


Appn.ove.cl  ^ok  pubtic  /(.eXcoAe;  dU&UbutZon  tmJUmittd. 


T158013 


Course-Keeping  With  Automatic  Control 

by 


Eduardo  Aguayo 
Lieutenant  Commander,    Peruvian  Navy 
B.S.,    Naval  Postgraduate  School,    1972 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGINEERING 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
Deceinber  1973 


A   - 


ABSTRACT 

A  ship  in  steering  is  considered  as  a  physical  mechanism  that  is 
forced  by  a  rudder  movement  to  produce  a  response.      Emphasis  is 
laid  upon  the  relation  between  the  forcing  and  the  response,    leaving 
aside  any  detailed  consideration  of  the  forces   concerned. 

Course-Keeping  with  Automatic  Control  techniques  is   studied 
following  this  concept.      Computer  programs  are  develoj^ed  to  simulate 
different  conditions.     Interpretation  of  the  results  is  made  to  evaluate 
the  different  methods  used. 
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I .     INTRODUCTION 

Problems  that  justify  the  use  of  Automatic  Control  for  ships: 

The  maneuverability  of  a  ship  is  determined  by  the  physical 
properties  of  the  ship  in  its  environment  and  the  physiological  and 
psychological  properties  of  its  navigator. 

A  successful  mission  of  a  vehicle  in  its  operating  environment  is 
a  function  of  its  handling  as  well  as  its  design.      The  relationship 
between  helmsmanship  and  the  steering  system  has  been  discussed  in 
detail  by  Abkowitz  [1],    Brard  [2],    Davidson  [3],    Nomoto  [4],    among 
others,    and  the  conclusions  indicate  that  the  hydrodynamic  design 
factors  alone  can  not  significantly  improve  the  handling  qualities  of 
ships . 

The  past  decade  will  be  known  as  the  decade  of  the  Super- Tankers 
Ships  were  made  larger  for  many  reasons,    one  of  which  was  to  im- 
prove operations  efficiency.      To  realize  the  full  benefit  of  the  larger 
size,    the  ship  had  to  be  steered  as   efficiently  as  possible.     Almost 
every  Super- Tanker  built  nowadays  is  course  unstable  and  the  control 
of  svich  a  system  has  been  of  great  concern.     We  can  see  easily  that 

i 

it  is  not  difficult  to  design  an  autopilot  if  enough  rate  control  is 
available  [5]. 

The  basis  for  the  need  of  Automatic   Control  was  a  minimum 
number  of  operator  controls,    minimum  rudder  orders  for  course- 
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keeping  and  automatic  course  correction  without  overshoot.      But  the 
main  reason  for  this  type  of  control  is  that  steering  an  unstable  ship 
by  manual  control  requires  a  helmsman  to  detect  a  very  small  amount 
of  rate  of  turn  and  to  react  to  this  in  time  with  the  proper  correction. 

Even  though  the  study  of  man  as  a  pilot  of  giant  tankers  is  rel- 
atively young,    yet  it  is   evident  that  any  pilot,    not  withstanding  his 
experience  and  seamanship,    has  some  physiological  and  psychological 
limitations   relevant  for  the  execution  of  his  task. 

Several  studies  carried  out  [6  and  7]  provided  us  with  some  in- 
sight.     First  a  purely  physiological  problem  is  a  limited  capacity  for 
the  perception  of  motion.     Very  slow  motions  cannot  be  perceived  by 
human  beings  as  is  the  case  with  the  hands  of  a  clock.      To  be  more 
specific,    the  smallest  yaw  velocity  human  beings   can  perceive  is 
about  1  minute  of  arc  per  second.      The  smallest  just  noticeable 
acceleration  or  deceleration  of  the  yaw  motion  occurs  when  the  yaw 
velocity  is  doubled  or  halved  within  5  seconds.     Analysis  of  ship  man- 
euverability experiments  on  full-scale  models   revealed  that  tankers 
above   100,  000  Ton.    move  so  slowly  d\iring  maneuvers  that  many  of 
the  accompanying  velocities  and  accelerations  are  not  perceptible 
for  the  man  on  the  bridge.      This  means  that  some  information  which 
is  indispensable  for  the  correct  execution  of  the  maneuvers  is  not 
available  for  the  navigator  due  to  his  own  physiological  limitations . 
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The  inertia  of  large  tankers  also  raises  a  problem  of  a  more 
psychological  nature.     It  is  well  known  to  pilots  and  Captains  that 
there  is  a  time  lag  between  a  rudder  command  and  the  reaction  of 
the  ship.      This  lag  can  be  easily  of  the  order  of  15  seconds.      Psy- 
chological research  on  tracking  behaviour  revealed  that  steering 
tasks   raise  difficulties  when  the  time  lag  is  longer  than  4  seconds. 
Also  related  to  inertia  is  the  ability  to  anticipate  motions  over  long 
time  intervals.      Even  very  simple  maneuvers  take  a  long  time  and 
should  therefore  be  initiated  a  long  time  beforehand.     It  is  however, 
very  difficult  to  anticipate  maneuvers  over  long  time  intervals  as 
human  beings  have  only  a  limited  capability  for  extrapolation  of  slow 
motions,    i.  e.  ,    the  navigator  can  hardly  predict  what  his  position 
will  be  in  10  minutes  and  therefore  is   not  able  to  initiate  the  neces- 
sary correction  at  the  required  moment. 

Still  another  factor  is  the  ability  to  discriminate  among  the 
different  components  of  the  ship's  motion,    those  due  to  wind,    current 
and  steering  movements  initiated  by  the  pilot. 

From  experiences  in  working  with  the  concept  of  the  time  con- 
stant of  a  ship,    conclusions   resulted  in  that  the  difficulty  in  steering 
depends  largely  on  this  parameter,    since  for  the  same  statistical 
characteristics,    the  unstable  ship  with  the  small  value  of  absolute 
time  constant,    changes  its  course  too  quickly  to  be  followed  by  the 
human  operator. 
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II.     COURSE- KEEPING 

A.      DEFINITION 

Course- Keeping  quality  is  defined  as  the  ability  of  a  ship  to  keep 
course  easily  on  a  seaway.      This  is  not  identical  to  stability  on  course, 
because  the  former  depends  upon  the  behaviour  of  the  steering  devices 
as  well  as  of  the  ship,    while  the  latter  is  a  thoroughly  passive  char- 
acter of  the  ship.     So  the  problem  of  course-keeping  should  be  treated 
by  analyzing  a  closed-loop  feedback  system  composed  of  a  ship,    an 
automatic  or  manual  steering  control  unit  and  steering  gear.     It  is 
possible,    however,    to  show  simply  that  good  stability  on  course 
yields  a  good  course- keeping  quality,    in  general,    in  this  form: 

1.  A  ship  that  is  more  stable  on  course  has  less  occasion  to 
deviate  significantly  from  its  course  because  the  stimulated  motion 
decays  more  quickly  even  with  the  rudder  amidship. 

2.  A  ship  that  is  more  stable  on  course  responds  more  quickly 
to  steering,    so  that  any  course  deviation  can  be  corrected  more 
easily. 

The  stability  on  course  is  related  to  Hie  decay  of  the  yaw  and 
sway  of  the  ship  with  a  rudder  amidship  after  being  disturbed  for  a 
short  while.     In  the  case  of  an  unstable  ship  this  results  in  a  negative 
damping  of  yaw  and  sway.     Stability  on  course  is  also  related  to  quick 
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response  to  steering.      This  is  of  particular  importance  in  considering 
the  course-keeping  quality  of  a  ship. 

B.      STEERING  QUALITIES  OF  SHIPS 

The  dynamic  stability  on  course  of  an  unsteered   ship  is  measured 
by  the  number  of  ship  lengths  traveled  by  a  stable  ship,    in  the  time 
required  to  reduce  an  accidental  deviation  from  undisturbed  motion 
in  a  straight  line  to   1/e  of  its  initial  value,    the  rudder  remaining  in 
the  amidship  position,    this  quantity  is  defined  as    1/lp    !  . 

A  negative  value  of  p.    moans  that  the  ship  is  dynamically  stable, 
i.  e.  ,    when  a  dynamically  stable  ship  moving  on  straight  course  is 
disturbed  slightly,    it  settles  down  on  a  new  straight  course  that  is 
close  to  the  original  one  (Figure   1).      The  greater  the  negative  magni- 
tude of  p.,    the  more  rapidly  the  ship  settles  on  its  new  course  and 
the  closer  the  new  course  is  to  the  original  one.     A  positive  value  of 
p,   means  that  the  ship  is  dynamically  unstable,    i.  e.  ,    when  a  dynam- 
ically unstable  ship  moving  on  straight  course  is  disturbed  slightly, 
it  will  go  into  a  steady  circuling  motion,    even  though  the  rudder  is 
held  amidship. 

The  primary  purpose  of  an  automatic  steering  device,    just  a.s  of 
a  human  helmsman,    is  to  compensate  for  disturbances,    making  the 
ship  maintain  a  prescribed  heading  and  thus  making  it  directionally 
stable. 
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FIGURE  1 


C.      UNSTABLE  HULL 

In  this  study  of  Automatic  Course-Keeping  we  will  be  working 
with  a  200,  000  DWT  super-tanker  of  the  following  characteristics: 

Length         310  meters 

Breadth         47.  16      meters 


Draft 


18.90      meters 


Steering  Quality  indices 


T        =     -269.3      seconds 


T 


K 


9. 3       seconds 

20.0      seconds 

-0. 0434    rad/sec 


Maximum  Rudder  deflection         30 


Maximum  Rudder  rate 


2.32  degrees /second 
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D.      AUTOMATIC  CONTROL  AND  COURSE. KEEPING 

One  of  the  functions  of  ship  control  is  to  maintain  a  ship's 
heading.     In  performing  this  function,    a  helmsman  deflects  the  rudder 
in  a  way  which  will  reduce  the  error  between  the  actual  and  desired 
heading  (Figure  2) 


v^  Actual 


FIGURE  2 

A  good  helmsman  will  not  only  deflect  the  rudder  in  response  to 
the  heading  error,    but  he  is  also  sensitive  to  the  angular  velocity  of 
the  ship  and  he  will  ease  off  the  rudder  and  apply  a  little  opposite 
rudder  in  order  to  prevent  overshooting  the  desired  heading.      It 
follows  that  an  automatic  control  should  also  be  responsive  to  control 
signals  measuring  both  error  and  angular  velocity.      Thus  a  rudder 
under  automatic  control  might  be  deflected  in  accordance  with  the 
following  linear  expression: 
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£2  angular  velocity  of  ship 
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The  basic  action  of  the  automatic  control  or   "autopilot"  is  called 
proportional  control  which  means  to  give  a  helm  angle  proportional 
to  the  amount  of  course  deviation  and  sometimes  also  proportional 
to  the  time  rate  of  the  deviation. 

In  general  the  complete  steering  system  can  be  represented  as 
in  figure  3 
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FIGURE  3 


The  response  of  a  ship  to  steering  or.    the  maneuverability  of  a 
ship  is  usually  described  by  a  set  of  equations  of  motion  of  side- 
drifting  velocity  and  turning-angular  rotation: 


jL 


1,,  4ft 


<U, 


a. 


ZZ 


r-     - 

% 

m 

■+ 

9 

h. 

/ 


It  is  however,    more  convenient  to  use  a  single  equation  of  motion 
describing  immediately  a  relation  between  turning  angular  motion  and 
steering,    because  control  signals  in  autopiloting  relates  only  to  turn- 
ing angular  motion. 
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111  •     TRANSFER  FUNCTION  IN  STEERING 

The  general  principle  of  linearizing  the  equation  of  motion  of  a 
ship  in  steering  may  be  summarized  as  follows: 

1.  The  coupling  of  surging  to  sway  and  yaw  is  ignored  because 
of  its  second- order  effect  on  the  latter  motions  [19]. 

2.  Hydrodynamic  forces  acting  upon  the  ship  are  expressed  as 
linear  functions  of  angle  of  drift,    i.  e.  ,    sway  velocity,    curvature  of 
the  ship's  path,    i.  e.  ,    yaw  angular  velocity,    and  angle  of  helm. 

3.  Hydro-inertial  forces  are  expressed  as  linear  functions  of 
the  acceleration  and  angular  acceleration  of  the  ship. 

Linear  equations  of  motion  of  this  kind  have  been  introduced  in 
a  number  of  articles  [8,  9,  10  and  17],    and  have  proved  promising  in 
interpreting  and  predicting  behaviours  of  ships  in  steering.      In  the 
notation  used  here,    they  are  written  as: 

(L/r)  (i»'  +  »i',)  ft  +y;  ?  -.  (L/r){*i'*m'K-Yr)e--Y!s 


in.  l 


J 


Using  these  equations  we  can  obtain  the  yawing  (  0    (t))  and  the 
swaying  (  p  (t))  in  response  to  any  movement  of  the  rudder  (  o    (t)). 

Now  limiting  our  interest  to  the  steering- to-yaw  response  which, 
in  general,    is  naturally  of  major  importance,    we  introduce  the 
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transfer  function  in  steering.      Taking  the  LaPlace  transform  of  both 
sides  of  the  equations  III-  1 ,    we  obtain: 

-  (  Vr)  («'+  M'y)  fir*)  +  C(L/r)  (*<'+  ^'y)5  +  Yji]  Pc*) 


-r/r)  (**'+**'*  -Yr)  &<*)  =    Y*    SC5) 

'^(z^+r^eco)  f[(v^x(2^9^rMjs  +(%/K]e<*) 

-   A/ ft     f*)      =    A/j     cTcsJ 


WJ 


'here  h    (0)  and    Q    (o)  denote      ft      and     0      respectively  at  t=0f 
namely  at  t   e  beginning  of  a  rudder  movement,    and  where 

fo     (5")     is  the  LaPlace  transform  of     (-> 
g    (S)     is  the  LaPlace  transform  of     (3 
<S>    (s )     is  the  LaPlace  transform  of     £ 
Taking  into  account  the  following  relation  at  t=0 

7 

(yirj  (Izz  +Jt*)Q\°)  +■(%) Aiy  9(o)  -a/^  fi(o)  -o 

we  can  eliminate    fi'(S)  and  3   [O)  from  the  transformed  equations  [19], 

to  yield  a  single  equation  that  describes  the  steering-to-yaw    re- 

1 
sponse,      that  is: 

G(S)-!LiJJtI*lL~  £cs)  -f  fi  Tz*,  +<rt  +  to] e(o)  v-  r,  r7_  0(o)  jn  2 


Since  later  in  this  study  we  are  going  to  be  dealing  with  auto- 
matic piloting,    and  the  control  signal  under  this  condition  relates 
only  to  turning  angular  motion. 
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where 


T, 


=  (V.) 


The  first  term  of  the  right  side  of  eq\iation  III- 2  corresponds  to  a 
motion  excited  by  steering  and  the  second  to  another  resulting  from 
the  initial  motion  at  t=0,    the  latter  disappears  if  the  ship  is  running 
straight  at  t  =  0  without  any  yaw  or  sway.      Thus  the  rational  function: 

B(s)  K  (  (+Tis) 

Jrs)     ~~    (  t+T,s)C^TiS) 

describes  the  response  behaviour  of  a  ship  in  terms  of  the  LaPlace 
transform.     We  may  call  it  the  transfer  function  of  the  ship  in  steering. 

Retransforming  equation  III-2,    we  obtain  a  single  differential 
equation  that  describes   ship  motion  just  as  does  the  original  equation 
of  motion  (Equations  III-  1 )  as  far  as  the  steering-to-yaw  response  is 
concerned,    that  is: 
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better  known  as  Nomoto's  Equation,    this  describes  directly  the  yaw 
response  of  a  ship.     The  four  coefficients  K,  T,,  T    ,    and  T„,    which 
are  composed  of  the  coefficients  of  the  original  equation  of  motion, 
constitute  a  set  of  characteristic  figures  representing  the  response 
behaviour  of  the  ship,    they  are  called  the  steering  quality  indices. 

A.      STABILITY  CRITERION  OF  AUTOPILOTING  BY  KOCHENBURGER'S 
METHOD 

A  control  loop  describing  autopiloting  of  a  ship  is  illustrated  in 
figvirc  3  in  the  form  of  a  block  diagram.     Since  the  stability  criterion 
of  Kochenburger  [11]  is  based  upon  whether  sinusoidal  signals  grow  or 
decay  in  circulating  through  the  loop,    it  is  necessary  first  to  obtain 
the  response  of  each  element  composing  the  system,    to  sinusoidal 
signals. 

The  response  of  a  ship  to  a  sinusoidal  signal  (that  is,    in  this    case, 
to  put  a  rudder  sinusoidally  to  both  sides  with  a  certain  frequencyfw) 
may  be  determined  through  the  equation  of  Nomoto: 

OCt)=  A(<*>)  So  ^  [>  fc  +  &  fWJ 

where 


Am  = 


__K   (  /   +  ceo  I ?_) 


(/  +  Cut  T,  )  c  /  ■(■  <-*a;  TO 
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and  where    $      is  the  amplitude  of  sinusoidal  steering.     A(oo»)  is  called 
an  amplitude  ratio  and   0    {uj)  a  phase  difference.      Both  of  them  are 
functions  of  frequency    w    only,    as  a  common  feature  of  linear  systems 

Next,    the  response  of  an  electrohydraulic  steering  gear,    which  is 
widely  used  for  most  present  ships,    may  be  described  by  the  following 


equation: 


where  T  is  a  time  constant  of  the  steering  gear. 

Then  we  get  a  description  of  the  response  of  a   steering  gear  to  a 
sinusoidal     ^    ,    as  follows: 


where 


£     ■-    A£  ceo;    X  g^u  £ 


U) 


&  «">] 


Ac  M  - 


/   +  t  co  T6 


(Pr-    M    =    An 


/     /    -+  c'cu7> 


Finally,    we  get  easily  the  response  character  of  a  compass  as 
follows,    because  it  may  be  considered  a  simple,  integrating  element 
transforming     0  into     ©    : 
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6>  =  Oo   *~  wt 

where 

Ac(y>)  - 


t  u> 


Oj 


£W  --   A?   r£    =    -  f 

Since  the  transmitting  character  through  all  linear  elements  has 
been  thus  obtained,    if  a  similar  character  of  the  reinaining  element, 
viz,    an  autopilot,    is  defined,    we  can  judge  whether  sinusoidal  signals 
grow  or  decay  in  circulating  through  the  control  loop.     It  is  impossible, 
however,    to  describe  the  response  of  an  autopilot  by  any  linear  differen- 
tial equation  and  then  to  obtain  its   response  character  in  the  foregoing 
manner,    if  considering  such  a  discontinuous  element  as  a  weather- 
adjust  mechanism,    which  is  used  to  modified  proportional  control  in 
order  to  avoid  frequent  steering  in  practical  application  for  rough 
seas. 

Fortunately,    however,    the  response  itself  of  an  autopilot  with  a 
weather-adjust  mechanism  to  a  sinusoidal  input  signal  (that  is,    in  this 
case  a  course  deviation)  may  be  obtained  easily  by  expanding  the  re- 
sponse into  a  Fourier  series  with  a  fundamental  frequency  that  is  the 
same  of  the  input  signal,    as  follows: 
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where      C  is  the  proportionality  constant  connecting  a 

*  course  deviation  to  a  helm  angle  to  be  called 

for, 


r£  is  the  amplitude  of  course  deviation,    and 

a    ,  a    ,a    ,    and  a         are  the  coefficients  of  Fourier  series  which 

1      2       V  4 

~*  may  be  obtained  through  the  usual  procedure 

of  Fourier  expansion  if  the  form  of    &       is  given. 
Considering  here  that  all  linear  elements  are  much  more  insensi- 
tive to  a  signal  with  the  higher  frequency,    we  can  neglect  all  the  higher 
frequency  terms.      This  is  the  approximation  of  Kochenburger,    and  its 
validity  depends  upon  how  much  the  higher  frequency  signal  decays 
through  the  linear  elements  and  also  how  much  higher  frequency  com- 
ponents  are  included  in  the  original     o    ;  viz,    how  much  the  original    <3 
resembles  a  pure  sinusoidal  form.     In  the  present  case,    Kochenburger 's 
approximation  may  be  fairly  valid  because  a  ship  is  q\iite  insensitive  to 
a  high-frequency  steering  because  of  her  large  inertia,    and  also  the 
basic  action  of  an  autopilot,    even  with  a  weather-adjust  mechanism,    is 
proportional  control  that  produces  a  sinusoidal  output  in  response  to 
a  sinusoidal  input. 

B.      FIRST  ATTEMPT  TO  STABILIZE  THE  SYSTEM 

The  autopilot,    the  steering  engine  and  the  ship  with  its   rudder  all 
form  different  components  of  a  closed  loop  system,    each  component 
characterized  by  its  transfer  function  or  the  complex  ratio  of  output  to 
input.      The  theory  for  such  control  systems  and  their  stability  has  been 
developed  in  electric  network  engineering,    and  it  is  natural  that  the 


25 


dynamic  problem  lends  itself  to  studies  in  analog  computers,    where 
each  component  is  represented  by  its   equivalent  electric  circuit.      The 
stability  of  the  closed  loop  system  may  be  judged  from  the  total  open 
loop  response  recorded  at  several  frequencies,    without  a  knowledge 
of  the  individual  transfer  functions. 

The  conditions  of  directional  stability  in  automatic  steering  along 
a  fixed  course  had  been  discussed  by  Minor  sky  [12]  by  means  of  the 
technique  of  added  derivatives,    applied  to  a  simplified  one-degree 
of  freedom  oscillation  and  including  several  types  of  position  and  rate 
control.     In  1946  Davidson  and  Schiff  [13]  took  a  large  step  towards  a 
better  understanding  of  the  interrelation  between  the  performance  of 
a  ship  on  a  straight  course  and  in  turning,    pointing  out  the  nonlinear  i- 
ties  in  the  behaviour  of  the  unstable  ship  and  including  an  interesting 
treatment  of  the  transients  when  entering  a  steady  turn.      They  also 
used  linear  theory  for  establishing  a  formula  for  the  radius  of  the 
turning  circle. 

The  work  of  Davidson  and  Schiff  did  much  to  stimulate  other 
authors.     Among  many  others,    a  paper  by  Williams  [14]  on  initial 
stage  motion  and  a  report  by  Schiff  and  Gimprich  [15],    who  studied  an 
automatic  control  system  where  the  rudder  angle  called  for  is  pro- 
portional to  a  combination  of  heading  deviation  and  rate  of  change  of 
heading,    and  which  has  a  behaviour  with  a  close  resemblance  to  the 
automatic  pilots  used  in  practice. 
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Most  modern  steering  engines  are  designed  to  move  the  rudder 
with  an  essentially  constant  speed,    the  rudder  turning  at  that  speed 
as  long  as  a  control  signal  is  transmitted  to  the  steering  engine.     A 
gyro  pilot  may  be  used  to  switch  on  this   signal,    calling  for  a  correct- 
ing rudder  at  a  certain  small  deviation  from  the  desired  heading, 
whereas  a  contact  on  the  rudder  may  stop  the  motion  at  a  suitable 
angle,    the  rudder  remains  in  this  position  until  the  ship  has  swung 
over  to  the  other  side  of  its   course  and  the  rudder  is  then  reversed. 
More  often  a  follow  up  mechanism  assures  a  "proportional  control" 
of  the  rudder.      Due  to  the  finite  rudder  speed,    these  systems  may  be 
self  exciting,    however,    and  excessive  oscillations  may  be  built  \ip. 

In  order  to  overcome  these  difficulties  the  automatic  pilot  must 
be  made  to  anticipate  the  motion  of  the  ship,    much  in  the  way  an 
experienced  helmsman  gives  an  auxiliary  rudder.      In  practice  this  is 
accomplished  by  means  of  a   feedback  of  rudder  motion  to  the  heading 
error  detector,    as  in  many  commercial  applications,    sometimes  also 
by  adding  some  kind  of  rate  of  change  of  heading  control,    similar  to 
the  pitch  rate  component  of  submarine  depth  control  systems.      In 
effect,    both  these  methods   correspond  to  a   "proportional  plus   rate 
control;"  the  first  one  often  incorporates  a.  non-linear  character  to 
the  system  [  16]. 

With  these  concepts  in  mind  we  are  going  to  try  to  stabilize  our 
unstable  system,    first  by  assuming  that  the  autopilot  is  just  a  gain  G, 
we  represent  the  block  diagram  in  figure  4 
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The  open  loop  transfer  function  for  this  svstem  is: 


Mrs)  .. 


(5    /<    C  /  -A    5  /  3 ) 


Sti  +  sr,)CusTz)  Cl+ste) 


From  experience  it  has  been  determined  that  T      ranges  from 
1  to  2,    and  a  good  choice  is    1.7  seconds,    so  the  only  parameter  avail- 
able to  make  the  system  stable  is  G. 

Since  there  is  one  pole  in  the  right  hand  plane,    given  by  T    ,    there 
is  no  sense  in  using  the  Bode  plot  to  determine  stability  values. 

1.       Root- Locus  Analysis 

Getting  the  characteristic  equation  and  arranging  it  in  the 
proper  form  we  obtain: 
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C£  :      1   +-AK(i< 


SO  +  sTjO  +  sTOO+su) 


* JLHI  te  -  j       =»     g     <-•<""*>  f"20*) =  -y 

Y(S)  j-  (/-lC<f.3s)(f+  9.3s)  (/  +  /.7S) 


Since  we  have  four  poles  there  are  four  separate  loci. 

There  are  three  asymptotes  due  to 

n  #  of  poles     =  4 

n-q  =  3 
q  /.'■  of  zeros    =1 

Center  of  asymptotes     -        _  _____I___~i_ — ■— -  —  —&.1I*) 

Angle  of  asymptotes  with  respect  to  the  real  axis: 

rr  +  v  2.TT  *$  ^  o,  i,  & t  ... .  n-9  -2 


71-9 

ir 


•P   =  a  C>o' 


0^2  /  VO* 

•$  ^  Z.  -Zoo" 

Froin  a  sketch  made  for  this  analysis  we  get  the  cross  point 

at  the  jw  axis  =  .  025.     So  G  at  the  crossing  point: 

^    '  I  S  j   t  S  f  o.  5-f  /  /  5  -  O.  oo*7 1  |  5  +  o.  lol  I 

Qy  >    _ . . ___.     3      Z 

Z  .oh  x  /c?"V  /  5f(5.o5( 

2.       Routh  Criterion 

To  find  a  better  set  of  numbers  we  refer  to  the  Routh 
criterion:  ■ 
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S      +  0.6,91    53    +O.OtS1^    (z.irr/o~"~  Z.ov  r/o*v£  )$  4-1. oi  x/S^C-o 
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where 


A     : 

6  = 


/  t    -  ■  8GP  Cr\ 


692 


A 


For  the  stability  limit  in  the  s     row,    B  must  be  equal  to 


zero  so: 


A  (-  '-«{%*/)  -<w  ('■<*  *>°*t)  = 


a 


after  manipulation: 


G    _   85   G  •!   240  =  0 


from  where 
82 


G  -     3        at  the  crossing  points 
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So  we  know  that  we  have  two  crossing  frequencies  in  the  jw 

2 
axis.     In  the  s      row 

2  -5 

As      +  1.02  x  10         G    =    0 

replacing  the  values  of  G,    and  after  all  the  mathematics  we  obtain: 


j     .0225 
iw  = 
J  j    .15 

values  that  agree  with  those  approximated  from  the  sketch. 

In  the  figure  5,    computer  output  of  the  subroutine  Root  Locus, 
we  obtain  the  open  loop  characteristic  of  this  system.     As  we  can  see 
the  system  can  be  considered  as  marginally  stable  due  primarily  to  the 
long  transient  that  can  be  predicted  from  the  root  locus  graph. 

Using  the  DSL  package  program  we  simulate  the  system, 
first  the  system  represented  by  the  block  diagram  of  figure  4,    and 
after  that  the  same  system  with  the  concept  of  "proportional  plus 
rate  control,  "  using  rate  of  turn  as  feedback,    as   represented  in 
figure  6 
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FIGURE  6 
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In  both  cases,    the  disturbance  was  simulated  as  a  step  of 
amplitude   1,    values  of  G  were  selected  between  12  and  36,    which  are 
the  ones  with  better  time  constant,    copy  of  both  computer  programs 
is  given  at  the  end  of  the  thesis. 

Figures  7-12  are  computer  outputs  of  the  first  condition, 
i.e.  ,    system  of  figure  4,    we  can  see  as   expected,    the  long  time 
required  to  settle  down  due  to  the  small  time  constant  of  the  G  values, 
but  even  so  we  can  see  in  figure  9  for  the  value  of  G  of  24.  2,    com- 
paratively a  short  time  for  reaching  steady  state. 

Figures    13-18  are  computer  outputs  of  the  second  condition, 
i.e.  ,    system  of  figure  6,    we  get  better  damped  results,    with  a  sig- 
nificant decrement  in  the  time  to  reach  steady  state  as  expected. 
Also  in  figure  15  we  can  see  that  comparatively  this  value  of  G  gives 
us  the  best  results.     So  we  can  state  that  the  value  of  G  =  24.  2  is  a 
good  value  to  stabilize  the  system  and  perhaps  the  best  under  this 
conditions. 
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I V .     BASIC  COURSE-KEEPING  REQUIREMENTS 

It  is  well  known  that  the  course-keeping  quality  and  the  turning 
quality  of  a  ship  are  contradictory  to  each  other,    and  that  it  is  im- 
possible to  improve  both  of  them  without  increasing  rudder  area  [4]. 

The  success  of  a  rudder,    designed  for  a  particular  ship,    is 
measured  by  the  degree  to  which  the  ship  achieves  a  desired  or  anti- 
cipated course-keeping  and  course- changing  ability.      Since  most  ships 
spend  most  of  their  operating  time  moving  ahead,    it  is  assumed  that 
maneuvering  performance  is   specified  only  for  the  ahead  direction. 

Course-Keeping  ability  may  be  specified  quantitatively  in  terms 
of  the  stability  index,    of  the  characteristics  of  the  Dieudonne  Spiral 

Maneuver  (see  Appendix  B),    or  of  the  range  of  rudder  angles  used 

2  r      . 

to  maintain  a  straight  course.        In  reference  [18],    it  is   suggested 

that  an  attempt  be  made  to  design  all  ships  for  a  stability  index  of 

zero  or  less,    but  it  is  recognized  that  this  may  not  be  practicable  for 

all  ships. 

A  major  reason  for  the  specification  of  a  level  of  control-fixed 

stability  for  all  ships  is  the  fact  that  excessive  controls -fixed 

instability  leads  to  excessive  use  of  the  rudder  to  maintain  a  straight 


2 

This  measure  of  Course-Keeping  ability  is  frequently  used  in 

maneuvering  in  restricted  waterways  or  in  rough  seas. 
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course.      This  in  turn  leads  to  increased  wear  on  the  rudder  system, 
decreased  ship  speed,    or  increased  fuel  consumption.      Controls- 
fixed  instability  also  leads  to  increased  difficulty  in  navigation  in 
restricted  waterways  and  in  following  seas  and  in  avoiding  collision 
with  passing  ships. 

A.      UNDER  STEADY  CONDITIONS 
1.       Calm  Water 

Among  the  factors  important  for  the  turning  motion  of  ships, 
there  are  transverse  forces  and  moments  acting  upon  ship's  body 
itself.      These  may  be  approximately  looked  upon  as  motion  in  infinite 
fluid  where  image  of  ship  is  considered  in  case  when  V  is  not  so  large 
and  the  effect  of  waves  by  the  ship  itself  with  Froxide  Number  below 


about  0.  2  or  0.  7     V     L/2d    [16]  is  negligible  or  when  motion  of  ship 
is  slow,    and  are  the  functions  of  V,    f->  ,     Q     .      This  will  be  sufficient 
for  turning  motion  of  ordinary  ship. 

The  disturbances  acting  tipon  the  system  can  be  classified 
according  to  their  influence  on  the  behaviour  of  the  system  as  follows: 

a.  Disturbances  that  cause  deviations  from  the  set  course. 

b.  Disturbances  which  affect  the  steering  characteristics  of 
the  ship. 

Wind  and  waves  belong  to  the  first  group,    to  the  second  class 
of  disturbances  belong  the  loading  of  the  ship,    the  depth  of  the  water, 
etc. 
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In  considering  the  calm  water  condition  we  are  assuming  no 
disturbance  at  all,    so,    under  these  conditions  we  expect  the  ship  to 
maintain  its  heading,    without  any  action  of  the  autopilot. 
2.       Steady  Wind 

Due  to  the  wind  pressure,    the  ship  will  be  acted  upon  by  a 
lift  and  drag  which  will  cause  the  ship  to  drift  at  a  drifting  angle.      By 
the  drifting,    the  hull  together  with  the  rudder  will  create  a  certain 
hydrodynamic  lift,    drag  and  moment  until  the  windward  component  of 
the  hydrodynamic  force  balances  with  the  wind  drag.     At  this   condition, 
the  total  summation  of  forces  including  ship's  thrust  is  zero.      How- 
ever,   the  summation  of  the  moment  including  moment  due  to  the  rudder 
at  maximum  helm  angle  will  not  necessarily  be  zero. 

If  the  summation  of  the  moment  is  not  zero,    and  is  in  a  direc- 
tion to  increase,    the  ship  will  have  a  tendency  to  turn  windward.     If 
the  moment  is  in  an  inverse  direction  the  ship  will  turn  leeward.     In 
both  cases,    the  ship  is   supposed  to  be  \mcontrollable.     If  the  summa- 
tion of  the  moment  is  zero  within  the  range  where  rudder  angle  does 
not  exceed  30  degrees,    the  ship  can  keep  its  heading  so  that  it  is 
supposed  to  be  controllable. 

B.      ACCURACY  NEED  FOR  COURSE-KEEPING 
1 .       Eoop  Gain 

The  most  important  parameter  to  be  considered  in  considering 
automatic  control  is  the  Loop  Gain,    and  its  value  depends  directly  on 
the  transient  response  of  the  system. 
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For  linear  systems,    the  most  commonly  used  correlation 

between  frequency  response  and  transient  response  is  that  correlation 

which  exists  between  the  height  of  the  resonance  peak  M        and  the 

pw 

height  of  the  peak  overshoot  of  the  step  response  M     for  a  second- 
order  system.     A  curve  relating  these  is  easily  calculated  [20]. 
Then,    for  a  known  second-order  frequency  response,    the  curve  is 
entered  with  Mpw  and  the  transient  peak  overshoot  M  *.  is  predicted 
exactly.     When  the  system  is  known  to  be  third-,    fourth-,    or  higher- 
order,    there  is  no  such  readily  available  correlation,    but  the  second- 
order  correlation  may  be  used  as  an  estimate.      This  gives  a  very 
accurate  estimate  if  the  higher-order  response  is  dominated  by  one 
pair  of  complex  roots,    a  less  accurate  estimate  if  dominance  is  not 
assured. 

We  know  that  in  order  to  make  a  system  stable  within  some 
degree  of  accuracy  we  need  to  have  complex  roots,    and  a  certain 
value  of  gain.      In  our  study  where  we  are  operating  with  a  giant  of 
200,  000  Tons.  ,    we  can  afford  to  have  a  certain  degree  of  inaccuracy, 
i.  e.  ,    1  to  2  degrees  of  heading  error,    without  consequences,    so  we 
can  go,    if  required,    to  have  real  roots  governing  our  system. 
2.       Loop  Type  Nximber 

Control  systems,    in  general,    are  required  to  have  certain 
operating  characteristics  which  are  determined  quantitatively  by  the 
specific  control  problem,    but  which  may  be  listed  qualitatively  the 
ones  we  are  interested  in  as  follows: 
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a.  The  system  must  be  acc\irate  in  steady  state. 

b.  The  system  must  be  stable. 

c.  The  system  must  regulate  against  disturbances. 

Each  problem  on  feedback  control  has  a  different  set  of  per- 
formance requirements  that  must  be  satisfied,    but  the  requirement 
that  is  common  to  virtually  all  problems  is  a  need  for  accuracy. 
What  is  meant  by  accuracy  depends  on  the  specific  physical  applica- 
tion,   and  no  convenient  broad  definition  is  available,    but  the  scope  of 
accuracy  requirements  is   readily  illustrated.      The  class  of  control 
system  we  are  interested  in  is  the  class  of  positioning  system  or 
servomechanisms.      These  sy steins  are  designed  to  change  the  output 
quantity  as  commanded  by  an  input  signal,    and  in  addition  are  required 
to  act  as  regulators  in  the  presence  of  output  disturbances.     The 
primary  consideration  is  usually  that  of  keeping  the  error  (difference 
between  commanded  output  and  actual  output)  less  than  a  specified 
amount  when  the  system  is  in  steady  state  with  no  load  disturbances. 
If  static  load  disturbances  are  anticipated  a  separate  specification 
for  steady  error  is  usually  given,    and  for  a  suddenly  applied  disturb- 
ance a  maximum  instantaneous   error  may  be  specified. 

For  positioning  systems  subject  to  command  input  signals, 
there  are  also  specifications  as  to  the  permissible  nature  of  the 
transient  response,    but  these  specifications  usually  refer  to  the 
permissible  time  duration  of  the  transient  and  to  the  permissible 
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nature  of  any  oscillations  during  the  transient  period.      They  seldom 
refer  directly  to  the  accuracy  of  the  system,    but  the  constraints 
which  they  place  on  the  system  design  may  make  it  very  difficult  to 
satisfy  the  accuracy  specifications. 

Systems   classified  as  positioning  systems  may  be  subjected 
to  various  types  of  input  commands.      These  may  be  well  defined  and 
representable  as   simple  mathematical  functions  or  they  may  be  rather 
complex  functions  of  time.      For  well  defined  inputs  the  system  usually 
(not  always)  reaches  a  finite  steady  state  condition  so  that  the  accu- 
racy of  the  system  is  readily  evaluated.      For  more  complex  inputs 
the  system  may  not  reach  a  clearly  defined  steady  state  condition, 
so  that  the  acc\iracy  must  be  evaluated  in  terms  of  bounds  and  a 
numerical  value  for   "steady  state  error"  may  not  exist. 

In  this   study  emphasis  is  placed  on  system  analysis  and 
design  for  cases  where  the  accuracy  problem  is  defined  in  terms  of 
deterministic  input  signals. 
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V .     DESIGN  OF  COMPE NSATION 

Having  been  able  to  obtain  the  best  value  for  the  gain  of  the 
autopilot  to  stabilize  the  unstable  model,    we  see  that  the  transient  is 
still  too  large  and  has  overshoots  (due  primarily  to  the  hydrodynamic 
coefficients  of  the  ship,    which  are  given,    as  we  may  recall  by  the 
steering  quality  indices). 

With  the  concepts  of  chapter  IV  we  are  going  to  approach  the 
design  of  the  compensation  using  the  Root- Locus  method  instead  of 
the  Bode  Plot,    since  having  a  pole  in  the  right  hand  plane  the  latter 
is  hard  to  use. 

A.      RATE  GYRO  COMPENSATOR 

Because  of  the  shape  of  the  Root  Locus  diagram  we  need  to 
improve  the  transient  response.      By  using  a  Rate  Gyro  we  introduce 
an  additional  zero  at  1/k,    and  this   relocates  the  roots  on  the  Root 
Locus  plot  at  a  different  gain  level.      Pvate  feedback  is  a  very  common 
means  of  increasing  equivalent  fiscous  damping  and  of  thus  improving 
system  transient  response.      Figures    19  and  20  are  a  block  diagram 
representation  of  heading  feedback,    the  original  circuit  and  the 
equivalent  respectively. 

The  new  characteristic  equation  is 
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Figures  2l_26  are  computer  outputs  of  the  subroutine  Root-Locus. 
They  give  the  different  system  behaviovir  for  the  values  of  k  selected. 
As  we  can  see  from  figures  24-26  in  order  to  have  a  good  time  constant 
we  need  a  large  value  of  G,    and  from,  figures  21_23  we  can  select 
very  small  values  of  G,    which  will  give  us   real  roots. 

B.      LEAD  TYPE  FILTER  COMPENSATOR 

Another  type  of  compensation  uses  a  filter  with  transfer  function 


of  the  form: 


S-f  £ 


<$  4  p 

having  high  pass  characteristics .     We  are  going  to  relocate,    also, 

the  roots  on  the  Root- Locus  diagram.      This  filter  is   just  a  cascaded 
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transfer  function  in  the  direct  path  between  input  and  output,    on  our 
system.     We  are  going  to  introduce  one  additional  pole  and  zero. 

Figure  27  shows  the  block  diagram  with  the  filter  in  the  system 


©« 


CF 


& 


-jPlLTC^ 


El- 


(■  t  !■  .'. 


FIGURE  27 


Skip 


-rtf^N^j-^ 


l*j 


J 


The  characteristic  equation  is: 

G  K  -k  Q-tsTz)  (s+i) 


I   + 


-    =o 


5  (  5  +  P)f  t+  s  T, J  (/  +571)  (f+S  jfcj 


by  selecting  the  values  of  z  and  p  as 


P 


.& 

.¥ 

,  oyy 

.</r 

.07J- 

J j 

.in 

U7 

.z 

2. 

,s* 

S. 

*/p 

7f 

,/ 

J.S 

./ 

X.  zz 

./ 

I.  3  3 

.  / 

.    f  6~5 

.  / 

.  6" 

.  I 

.  1 

we  obtain  the  5th  order  polynomials  to  use  in  the  Root- Locus  plot 
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Figures  28-33  are  computer  outputs  of  the  subroutine  Root-Locus. 
They  give  the  different  system  behaviour  for  the  values  of  z  and  p 
selected.     We  can  see  from  figures  31-33  that  in  order  to  have  a 
good  time  constant  we  need  a  large  value  of  G,    while  from  figures 
28-30  we  can  select  very  small  values  of  G,    which  will  give  us   real 
roots. 

Having  obtained  the  computer  outputs  for  both  types  of  compensa- 
tion,   we  can  observe  a  similarity  in  the  system  behaviour  for  same 
root  location.     So  we  can  predict  the  same  transient  for  step  inputs. 
However,    we  need  to  find  an  optimum  value  of  k  in  the  case  of  the 
Rate  Gyro  and  of  z  for  the  Filter,    in  order  to  expand  the  little  loop 
in  figures  2l_23  and  28-30.      This  optimum  is  the  best  value  possible 
to  obtain  and  it  is  a  compromise  between  a  high  gain  and  a  good  rudder 
operation,    bearing  in  mind  that  for  this  to  occur  we  must  consider 
maximum  deflection  and  rate  for  the  rudder  of  this   ship.     Since  for  a 
quick  response  we  need  a  high  gain,    but  this  could  demand  an  exces- 
sive operation  of  the  rudder  and  this  could  be  bad  for  the  following 
reasons : 

1.  To  get  quick  response  we  need  faster  rudder  angle  and  a 
greater  rudder  angle. 

2.  Faster  responses  may  cause  accelerations  that  are  too 
sharp  for  personnel  on  ship. 
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So  the  next  step  will  be  to  try  to  find  the  optimum  value  men- 
tioned above.      Right  off  hand  we  can  consider  the  system  behaviours 
of  fig\ires  25,    26,    32  and  33,    are  very  much  similar  to  our  original 
system  without  compensation,    so  we  can  disregard  those  four  systems, 
and  now  we  have  four  different  systems  for  the  Rate  Gyro  and  Filter 
compensators   respectively. 

We  are  going  to  simulate  a,  disturbance  that  produces  a  rate  of 
turn  of  .2  degrees /second,    a  value  that  has  been  found  by  experience 
can  easily  be  reproduced  on  this  tanker,    requiring  a  correcting  rudder 
angle  of  10  to  20  degrees. 

With  this  external  disturbance  represented  as  a  step,    we  examined 
our  eight  compensated  systems.      By  trial  and  error  method  we  check      U^jL^^ 


for  the  optimum  value  required  and  we  end  up  with  a  1/k  -  .  041  for 
the  Rate  Gyro,    and  a  z  =  .041  for  the   Filter,    which  in  essence  are 
the  same  values  we  have  for  the  case  of  1/k  =  .  04  and  z  =  .  04  for  the 
Rate  Gyro  and  Filter  respectively,    so  now  remain  these  two  systems 
to  be  studied  for  the  best  gain.      In  all  the  eight  initial  cases  we  ob- 
served a  common  characteristic,    with  minor  variations,    that  for  a 
small  variation  of  heading,    a  sudden  change  in  rudder  occurs  which 
causes  a  very  quick  response  and  we  may  suspect  an  excess  in  both 
rudder  parameters. 

Figures  34-37  are  computer  outputs  for  the  system  compensated 
with  the  filter,    we  have  chosen  values  of  gain  from  the  root-locus 
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graph,    for  four  different  damping  coefficients,    9.42,    15.87,    23.67 
and  3.  5   respectively.     We  can  observe  from  the  first  three  figures 
that  a  quick  response  is  obtained  and  steady  state  is  obtained  in  a 
very  short  time,    but  this  coiild  exceed  the  restrictions  of  the  rudder. 

In  figure  37  we  see  a  moderate  response  and  we  do  not  reach  the 
steady  state  at  the  end  of  the  simulated  time,    which  is  only  5  minutes, 
as   expected  for  a  ship  of  characteristics   such  as  this.     We  need  to 
verify  if  there  is  any  violation  of  the  rudder  restrictions. 

With  the  help  of  computer  simulation,    a  program  in  DSL,    which 
is  given  at  the  end  of  the  thesis,    is  used  to  simulate  limiters  for  both 
parameters.      In  order  to  introduce  the  limiting  value  of  Rudder  Rate, 
we  need  to  replace,    in  the  block  diagram,    the  steering  Gear  transfer 
function  for  a  circuit  equivalent,    given  in  figure  3  8: 
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where 


K     =  1/T 

3  E 


With  the  program  above  mentioned  and  the  values  of  gain  selected 
before  we  obtain  the  figures  39-42,    computer  outputs,    where  we  see 
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in  the  first  figures  that  even  if  the  system  is   stabilized,    there  is  a 
portion  of  the  response  where  we  exceed  the  rudder  deflection,    while 
in  figure  42  we  see  a  smooth  settle  down  without  any   clipping.     Will 
be  a  good  idea  to  try  values  of  gain  within  this   small  range  of  gain 
values. 

Figures  43-48  are  computer  outputs  for  the  system  above  men- 
tioned with  values  of  gain  of  6.  ,    5.3,    4.  6,    3.5,    2.5  and  2.     We  can 
observe  in  figure  44  for  a  gain  value  of  5.  3  we  obtain  a  good  settling 
time  with  the  smaller  steady  state  error  of  only  .  8  of  degree  out  of 
course  which  is   excellent,    but  we  have  some  rudder  action  even 
after  settling  time.     In  figure  45  for  a  gain  value  of  4.  6  we  obtain  a 
steady  state  error  of  1.  03  degree  out  of  course,    with  a  settling  time 
of  about  2  more  minutes  than  for  the  former  case,    but  we  get  smaller 
rudder  deflections.     So  we  can  conclude  that  this  value  of  4.  6  will  be 
the  best  value  for  this   system. 

Now  considering  the  Rate  Gyro  having  a  zero  at  1/k  =  .  04,    with 
the  same  disturbance,    we  obtain  the  computer  outputs  of  figures  49-52, 
for  values  of  gain  selected  from  the  Root-Locus   graph,    for  four  dif- 
ferent damping  coefficient,    22.26,    67.28,    35.29  and  3.5   respectively, 
we  observe  almost  the  same  results  as  in  the  case  of  the   Filter,    i.  e.  , 
the  first  three  figures  evidently  violate  the  rudder  restrictions  and 
only  in  figure  51  for  a  gain  value  of  3.  5  gives  us  a  good  system 
operating  under  real  conditions.      Figures  53-58  are  computer  outputs 
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for  this  system,    for  gain  values  of  6.  2,    5.  6,    5.  ,    4.  4,    3.  8,    and  3.  2, 
following  the  same  considerations  as  with  the  Filter.     We  can  see 
from  these  figures  that  in  all  cases  the  system  has  increased  its 
transient,    being  the  best  value,    6.2  as  observed  in  figure  53.     In  the 
event  of  increasing  the  gain  to  6.  7  we  start  obtaining  clipping  in  the 
rudder  action,    so  this  value  of  gain  is  chosen  as  the  best  suited  for 
this  type  of  system. 

From  comparison  between  these  two  systems  we  can  select  the 
one  with  Filter  compensation,    first,    because  of  the  best  final  results, 
i.  e.  ,    less  amount  of  rudder  required,    better  settling  time,    and 
secondly  for  the  reliability  of  the  system,    since  the  Rate  Gyro  could 
not  be  suitable  for  marine  use  due  to  this  overhaul  interval. 
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VI.     CONCLUSIONS 

A.      RESULTS 

An  optimal  value  of  gain  can  always  be  obtained  using  the  Root- 
Locus  method  for  any  given  plant.      Both  negative  and  positive  feed- 
back can  be  used  to  achieve  the  desired  results  but  only  negative 
position  feedback  should  be  used  to  obtain  a  stable  system. 

With  the  concepts  of  Chapter  IV,    the  requirements  for  a  basic 
Course-Keeping,    and  our  best  value  for  the  system  selected  as  the 
best,    we  verify  with  a  low  order  disturbance  represented  by  an  initial 
condition  of   O    »    that  the  actual  heading  in  the  steady  state  is  equal 
to  the  reference  heading  or  which  is  the  same  the  steady  state  error 
is   equal  to  zero.     We  already  have  observed  that  for  a  load  disturbance 
represented  by  the  step  input  we  obtain  a  small  steady  state  error. 
So  we  have  reasons  to  believe  that  we  are  working  with  a  type  zero 
system  and,    if  it  is  so,    do  we  need  to  change  the  loop  type  number? 

By  doing  an  analysis  of  the  feedback  loop  we  can  easily  verify 
that  indeed  we  have  a  type  zero  system,    and  using  the  final  value 
theorem  for  this  linear  system,    the  evaluation  of  steady  state  error 
for  the  deterministic  input  is  easily  accomplished.      In  this  case 
where  steady  state  error  is  finite  its  magnitude  is  determined  by  the 
reciprocal  of  the  gain  transfer  function  for  the  direct  path  between 
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input,  and  output,    this  magnitude  is  called  an  error  constant,    also 
called  the  DC  gain  or  zero  frequency  gain,    because  it  is  the  numerical 
value  obtained  for  the  transfer  function  by  deleting  the  s  factors  and 
substituting  s  =  0.      Following  this   reasoning  and  with  the  final  results 
obtained  we  see  that  this  meets  our  requirements  for  the  stable 
operation  of  the  system  with  the  degree  of  accuracy  that  can  be  con- 
sidered as  satisfactory. 

The  purpose  of  this  thesis,    which  was  to  study  the  basis  of  the 
Course-Keeping  with  Automatic  Control,    is   reached  at  this  stage. 
It  is  good  to  point  out  that  further  studies   could  be  made  by   studying 
the  sea  effects,    how  automatic  steering  and  the  yawing  of  ships  is 
affected  in  rough  seas. 
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APPENDIX  A 
LIST  OF  SYMBOLS 

d  draught  of  a  ship,  mean  of  fore  and  aft  if  no  remark. 

I  moment  of  inertia  of  a  ship  about  a  vertical  axis  through 

her  centre  of  gravity. 


x 


x 


N 


NJ, 


N', 


Ni 


■/■£■  *<■ 


J  moment  of  inertia  of  additional  mass  about  a  vertical  axis 

zz 

through  the  centre  of  lateral  additional  mass. 


L  length  of  a  ship,    between  perpendiculars,    if  no  remark. 

m  mass  of  a  ship. 

m  longitudinal  additional  mass. 


™X/L 


/£& 


m  lateral  additional  mass. 

y 

N  hydrodynamic  moment  abotit  a  vertical  axis  through  the  C.  G. 

of  a  ship,    positive  for  standard  turning  moment. 


yaw  angular  velocity,    positive  to  starboard,    identical  with  0 
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r'  r/(V/L)    or    L/R 

R  turning  radius,    steady  or  instantaneous 

V  ship  speed  in  meters /second 

Y  lateral  component  of  hydrodynamic  force,    positive  to  star, 
board. 

P 


y^     ^/7b/3 


/3  angle  of  drift,    positive  to  port. 

£  angle  of  helm,    positive  to  starboard. 

Q  turning  angular  velocity,    positive  to  starboard,    identical  to  r, 
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APPENDIX  B 
THE  DIEUDONNE  SPIRAL  MANEUVER 

The  spiral  maneuver  consists  of  the  following: 

1.  The  ship  is   "steadied"  on  a  straight  course  at  a  preselected 
speed  and  held  on  this  course  and  speed  for  about  1  minute.     Once  a 
steady  speed  is  established,    the  power  plant  controls  are  not  manip- 
ulated for  the  duration  of  the  maneuver. 

2.  After  about  1  minute,    the  rudder  is  turned  to  an  angle,    £     , 
of  about  15  degrees  and  held  until  the  rate  of  change  of  yaw  angle 
maintains  a  constant  value  for  about  1  minute. 

3.  The  rudder  angle  is  then  decreased  by  a  small  amount 
(about  5  degrees)  and  held  fixed  again  until  a  new  value  of    (3     is 
achieved  and  is  constant  for  several  ininutes. 

4.  The  foregoing  procedure  is  repeated  for  different  rudder 
angles  changed  by  small  increments  from,    say,    large  starboard     • 
values  to  large  port  values  and  back  again  to  large  starboard  values. 

The  numerical  measures  obtained  from  the  preceding  spiral 
maneuver  are  the  steady  yawing  rates  as  a  function  of  rudder  angle. 
A  plot  of  these  values  is  indicative  of  the  stability  characteristics 
of  a  ship.      For  example,    if  the  plot  is  a  single  line  going  from  star- 
board rudder  to  port  and  back  again,    as   shown  for  ship  A  in  figure   1-B, 
the  ship  possesses  controls-fixed,    straight-line  stability;  that  is,    it 
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has  a  negative  stability  index.     If,    however,    the  plot  consists  of  two 
branches  joined  together  to  form  a   "hysteresis"  loop,    as   shown  in 
ship  B  of  figure   1_B,    the  ship  is  unstable;  that  is,    it  has  a  positive 
stability  index.     In  addition,    the  height  and  width  of  the  loop  are 
numerical  measures  of  the  degree  of  instability;  the  larger  the  loop 
the  more  unstable  the  ship.     The  slope  of  the  yaw- rate  curve  at  zero 
rudder  angle  is  a  measure  of  the  degree  of  stability  or  instability. 
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